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ABSTRACT 

ANALYSIS  OF  LASER  INDUCED  SKIN  BURNS  BY 
A  DAMAGE  INTEGRAL  MODEL 


OBJECTIVE 


To  determine  the  applicability  cf  a  proposed  mathematical  model 
for  thermal  injury  predictions  to  experimental  data  from  a  study  of 
CO2  laser  radiation  induced  skin  bums. 

METHODS 


The  exposure  values  for  the  production  of  uniform  full  thickness 
epithelial  burns  were  used  to  determine  the  inactivation  rate  constants 
which  best  fit  the  data  over  the  exposure  time  range  0.2  to  5  sec.  The 
appropriate  time- temperature  values  were  calculated  by  heat  flow  equa¬ 
tions  based  on  an  assumed  isotropic  semi-infinite  opaque  solid. 

CONCLUSIONS 


The  data  within  the  exposure  time  range  of  0.2  to  5  sec  are  accu¬ 
rately  described  by  the  model  and  can  be  quantitatively  predicted  by  the 
equation 
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64 
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where  is  the  temperature  at  time  t  and  0  is  the  damage  integral  with 
an  arbitrarily  assigned  value  of  1.0.  For  times  in  excess  of  5  sec  the 
values  predicted  by  the  model  were  higher  than  those  determined  experi¬ 
mentally. 
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ANALYSIS  OF  LASER  INDUCED  SKIN  BURNS  BY 
A  DAMAGE  INTEGRAL  MODEL 


INTRODUCTION 


Programs  set  up  to  evaluate  the  potential  hazards  of  lasers  involv¬ 
ed  in  military  operations  have  been  largely  empirical  in  nature.  For  a 
given  frequency  of  radiation,  the  appropriate  combinations  of  power 
density  and  exposure  time  required  to  produce  some  designated  level  of 
injury  have  been  determined  by  trial  and  error.  When  the  large  number 
of  possible  combinations  of  wavelength,  power  density,  exposure  time, 
and  the  physical,  chemical,  and  biological  parameters  of  target  tissue 
are  considered,  the  need  for  a  more  basic  approach  becomes  quite  evident. 
Such  an  approach  would  limit,  to  a  reasonable  dimension,  the  number  of 
experiments  that  are  required  to  define  the  hazards  and  safety  limits 
associated  with  laser  use.  An  understanding  of  the  underlying  physical 
and  biological  mechanisms  involved  in  injuries  induced  by  laser  radia¬ 
tion  would  permit  the  extrapolation  of  data  from  a  limited  number  of 
experiments  to  other  appropriate  situations. 

It  is  generally  conceded  that  injuries  resulting  from  the  use  of 
lasers  in  military  applications  will  be  thermal  in  nature;  that  is,  they 
will  arise  from  temperature  increases  in  the  exposed  tissue  resulting 
from  the  degradation  of  the  radiant  laser  energy  to  thermal  energy.  The 
basic  principles  involved  in  the  analysis  of  thermally  induced  injury  in 
cutaneous  tissue  were  defined  in  a  series  of  papers  by  Henriques  and 
Moritz.  Henriques  (1)  considered  the  development  of  thermal  injuries 
to  be  time  dependent,  a  rate  process  that  is  common  to  most  temperature- 
sensitive  physical  or  chemical  systems.  He  suggested  that  by  means  of  a 
general  expression  applicable  to  all  rate  processes,  as  used  in  chemical 
kinetic  theory,  a  mathematical  equation  could  be  developed  to  predict 
quantitatively  the  time- temperature  thresholds  of  tissue  injury.  By 
integrating  the  injury  rate  process  over  the  time  period  involved,  the 
contributions  from  each  time  increment  would  be  summed  to  yield  the 
total  damage.  The  proposal  has  since  become  known  as  the  damage  inte¬ 
gral  model.  In  order  to  calculate  accurately  the  total  accrued  damage, 
the  precise  time- temperature  history  imposed  upon  and  the  thermal  sen¬ 
sitivity  of  the  biological  system  under  consideration  must  be  known. 
Recently  published  data  by  Brownell  et  al  (2)  for  CO2  laser  induced 
porcine  cutaneous  burns  offer  the  opportunity  to  apply  the  idea  of  the 
damage  integral  model  to  injuries  induced  by  exposure  to  laser  radia¬ 
tion  in  the  time  range  0.2  to  20  sec.  The  model  used  is  essentially 
that  proposed  by  Henriques  (1)  and  subsequently  expanded  by  Fugitt  (3) . 

It  is  not  the  intent  of  this  report  to  define  accurately  either  the 
time- temperature  relations  in  the  tissue  01  the  damage  rate  function 
involved,  but  only  to  determine  the  compatibility  of  the  proposed  mod¬ 
el  with  the  experimental  data. 

DAMAGE  INTEGRAL  MODEL 
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It  is  assumed  that  the  kinetics  of  the  thermal  inactivation  of  the 
tissue  cells  or  protein  component  of  the  tissue  that  is  used  as  an  in¬ 
dicator  of  thermal  damage  can  be  regarded  as  a  first  order  reaction  at 
all  temperatures  under  consideration.  The  rate  at  which  the  reaction 
will  proceed  can  then  be  given. 


dn 

at  = 


-  k'n 


(1) 


where  dn  is  the  number  of  elements  inactivated  in  the  time  increment 
dt,  k'  is  the  specific  reaction  rate  constant  for  a  given  temperature 
and  environmental  state  and  n  is  the  number  of  elements  at  any  instant. 

Upon  rearrangement  the  integral  of  equation  (1)  is 


or 
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where  n  represents  the  number  of  active  elements  remaining  at  time  t 
and  n0  represents  the  number  of  elements  at  time  t0.  The  specific  re¬ 
action  rate  k'  can  be  regarded  as  a  constant  only  for  a  given  tempera¬ 
ture  and  it  will  vary  as  the  temperature  varies.  If  equation  (2)  is 
integrated  over  the  time  interval  t0  to  t,  the  resulting  equation  will 
yield  a  numerical  value  for  the  surviving  fraction  for  the  temperature 
and  time  of  interest  by  utilizing  the  appropriate  k'. 

During  the  interval  in  which  a  thermal  injury  is  induced  by  radia¬ 
tion  (a  type  such  as  is  under  consideration)  the  temperature  is  not 
constant  but  varies  continually  during  the  episode.  Following  the 
onset  of  exposure,  the  tissue  temperature  rises  until  the  end  or  near 
the  end  of  the  exposure  period,  subsequently  it  falls  to  the  initial 
temperature.  Therefore,  to  complete  the  mathematical  description  of 
the  process  an  expression  is  needed  to  relate  the  inactivation  rate 
constant  to  temperature. 
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If  the  reaction  leading  to  thermal  injury  is  similar  to  that  of 
many  chemical  and  physical  rate  processes,  the  thermal  inactivation 
rate  constant  should  vary  with  temperature  as 

k,  .  ta-»/(T  •  2731  (4, 

where  A  and  B  are  constants.  If  the  reaction  meets  the  necessary  therm¬ 
odynamic  criteria,  B  is  equivalent  to  A  E/R  where  AE  is  the  energy  of 
activation  and  R  is  the  gas  constant;  the  equation  becomes  equivalent 
to  the  integrated  form  of  the  Arrhenius  equation.  In  this  case  A  and  B 
will  be  considered  experimentally  determinable  constants  with  no  special 
implications. 

Combining  equations  (3)  and  (4)  and  rearranging  yields 


in  Cn0/n)  =  A  j‘  *  2,3>  dt  (3) 

to 

where  Tt  is  the  temperature  at  time  t. 

It  is  practically  impossible  to  determine  the  number  of  cells  or 
protein  molecules  inactivated  in  the  skin  following  a  thermal  episode. 
Therefore,  an  arbitrary  level  of  tissue  injury  as  judged  by  histological 
and  histochemical  techniques  must  replace  the  quantitative  measure  as 
specified  by  the  express;on  In  (n0/n) .  Equation  (5)  then  becomes,  as 
suggested  by  Henriques  , 

<1  .  A  f  *  27S>  dt  (6) 


where  Q  is  an  arbitrarily  chosen  level  of  tissue  injury  and  is  desig¬ 
nated  as  the  damage  integral.  Equation  (5)  allows  calculation  of  the 
sm  of  the  damage  over  a  specified  time  interval  and  temperature.  A 
plot  of  Henriques’  data  (1),  where  he  determined  the  shortest  time  at 
which  a  constant  predetermined  tempeiature  produced  transepi dermal  ne¬ 
crosis,  demonstrates  a  good  fit  to  equation  (5)  for  conditions  under 
which  the  temperature  at  the  basal  layer  of  the  epidermis  could  reason¬ 
ably  be  assumed  to  be  constant. 

In  order  to  use  equation  (5)  to  sum  the  damage  accumulated  during 
a  thermal  episode  in  which  the  temperature  is  varying,  it  is  necessary 
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to  define  a  function  relating  temperature  to  the  time  of  the  thermal 
episode.  The  opaque  semi-infinite  model,  described  bv  Carslaw  and 
Jaeger  (7)  and  extended  by  Davis  (8,  9)  to  apply  to  problems  involving 
radiant  heating  of  tissues,  is  considered  the  most  appropriate  to  be 
used  in  the  calculation  of  time- temperature  relationships  at  variou.- 
tissue  depths. 

OPAQUE  SEMI- INFINITE  THERMAL  MODEL 

The  data  to  be  analyzed  is  taken  from  a  study  using  CO2  radiation 
with  a  wavelength  of  10.6  v  as  the  source  of  heat.  It  is  assuned  that 
the  porcine  skin  is  opaque  to  radiant  energy  of  this  wavelength  and  the 
radiation  is  absorbed  at  the  surface  of  the  skin.  The  half- layer  value 
of  skin  and  other  wet  tissue  for  10.6  u  radiation  is  approximately  10*3 
cm  with  an  absorption  coefficient  of  about  7C0  reciprocal  cm  (5,  6). 

The  following  tssuuptions  are  pertinent  to  the  opaque  semi-infinite 
model  as  used:  the  skin  is  a  serai-infinite  isotropic  receiver,  initially 
at  a  uniform  temperature  throughout;  the  thermal  properties  of  the  skin 
are  constant  and  do  not  change  with  temperature;  the  surface  is  perfect¬ 
ly  insulated  and  there  are  no  radiation  losses;  the  radiation  input  is 
a  square  wave,  normal  to  and  uniform  over  the  surface  of  the  tissue; 
heat  flow  is  unidirectional,  being  perpendicular  to  the  skin  surface; 
and  the  tissue  is  opaque  to  the  radiation,  which  is  absorbed  at  the 
surface. 

Under  these  conditions,  the  equations  describing  the  temperature 
change,  as  given  by  Davis  (8,  9),  are  as  follows 
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I  =  temperature  rise  (*C) 
.  =  depth  (cm) 
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t  =  time  (sec) 

- 1  -2 

I  =  irradiance,  absorbed  (joules  sec  cm  ) 
n  ■  exposure  time  (sec) 

2  -4  -2  -1 

jj  =  thermal  inertia,  kpc  (cal  cm  dig  cm  ) 

a  =  thermal  diffusivity,  k/pc  (cm^  sec  *) 

k  =  thermal  conductivity  (cal  sec  *  cm  *  dig  *) 

c  =  specific  heat  capacity  (cal  gm  *  dig  *) 

p  =  density  (gm  cm  "*) 

erfc  =  complementary  error  function 

The  temperature  change  during  the  exposure  period  is  given  by  equation 
(7)  and  following  the  termination  of  radiation  exposure  by  equation  (8). 
The  numerical  values  of  the  thermal  constants  used  in  these  calculations, 
u  =  11.8  x  10"4  and  a  =  8.6  x  10-4,  were  derived  experimentally  by  Davis 
(8,  9)  for  porcine  skin.  A  temperature  of  36°C  was  assumed  for  the  tem¬ 
perature  at  time  zero  in  all  cases. 

The  temperatures,  calculated  as  a  function  of  time,  were  then  used 
in  equation  (6)  to  determine  the  numerical  value  of  the  damage  integral 
for  each  exposure  time.  Sixty  to  80  equal  time  intervals  were  summed 
over  the  time  during  which  tha  damage  rate  was  greater  than  1/100  of 
the  maximum  damage  rate  for  each  case.  Damage  integrals  were  deter¬ 
mined  only  for  those  cases  where  the  exposure  time  was  less  than  in  sec. 
It  was  thought  that  the  dermal  capillary  bed  may  act  as  a  heat  sink  for 
longer  periods  of  time,  thereby  invalidating  the  assumption  of  the  semi¬ 
infinite  isotropic  feature  of  the  thermal  model.  Constants  A  and  Bin 
equation  (6)  were  varied  to  give  the  closest  fit  to  a  damage  integral 
for  all  cases  within  the  exposure  time  range  of  .2  to  10  sec. 

RESULTS  AND  DISCUSSION 

The  data  obtained  from  the  gross  evaluation  of  2,288  porcine  skin 
burns  induced  by  exposure  to  CO2  laser  radiation  (2)  is  presented  graph¬ 
ically  in  Figure  1  (next  page).  The  semi logari thmic  plot  gives  the 
radiant  energy  density  in  jouies/cm-  that  produces  at  least  a  minimally 
detectable  erythema  50^  of  the  time.  This  was  uesignated  as  a  1-2  bum. 

Selected  biopsies  were  made  in  order  to  establish  a  relationship 
between’ the  surface  appearance  of  the  burn  and  the  depth  of  damage  as 
judged  histologically  and  hisf  ocheraical ly  (3).  Of  the  1-2  burns  ex¬ 
amined  microscopically,  approximately  ~0*  weie  uniform  full  thickness 
epithelial  lesions,  7\  included  dermal  damage,  3nd  the  remainder  had 


5 


only  partial  epithelial  damage.  Although  the  data  suggest  that  the 
depth  of  damage  may  increase  slightly  as  the  power  density  increases, 
they  were  not  sufficient  to  establish  this  effect.  The  average  thick¬ 
ness  of  the  epithelium  was  75  u. 

Based  on  the  biopsies,  a  1-2  burn  is  considered  equivalent  to  a 
histologically  defined  uniform  full  thickness  epithelial  lesion.  There¬ 
fore,  the  data  in  Figure  1  represent  the  radiant  energy  density  values 
with  a  501  probability  of  producing  thermal  damage  to  a  depth  of  0.075 
cm  in  porcine  skin.  This  is  designated  as  the  damage  level. 
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Fig.  1.  The  incident  energy  density  required  to 
produce  visually  detectable  erythema  50%  of  the 
time  as  a  function  of  the  exposure  time. 


By  fixing  the  thermal  constants,  a  and  y,  and  arbitrarily  assigning 
a  numerical  value  of  1.0  to  the  damage  integral  SI,  it  is  possible  to 
find  the  constants  A  and  B  for  the  damage  rate  function  that  best  fit 
all  the  data  from  .2  to  10  sec  exposure  time.  The  following  values  were 
found  to  give  a  reasonably  constant  damage  integral  of  1.0  for  nearly 
all  the  cases  examined: 
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In  A  =  147.518 


B  =  49,110 

The  total  damage  integral,  as  calculated  for  each  exposure  time- 
power  density  combination  is  given  in  column  6,  Table  1.  With  the  ex¬ 
ception  of  exposure  times  .37,  .71,  and  7.4  sec,  the  calculated  damage 
integrals  are  quite  close  to  the  assigned  value  of  1.0. 


TABLE  1 

Temperature  and  Damage  Integral  Values  Derived  From  Experimental  Data 


Exposure 

Time 

(Sec) 

Irradiance 

(watts/cm^) 

Max. 

Temp. 

(°C) 

Time  of 

Max.  Temp. 
(Sec) 

Damage 
Integral 
to  end  of 
Exposure 

Total 

Damage 

Integral 

.22 

13.6 

56.2 

.230 

0.55 

1.05 

.30 

10.6 

65 . 6 

.309 

0.51 

0.98 

.37 

7.6 

60.6 

.378 

0.08 

0.15 

9.0* 

65.2 

.378 

0.50 

0.94 

.71 

4.7 

59.6 

.717 

0.09 

0.16 

5.6* 

64.1 

.717 

0.60 

0.99 

1.3 

3.7 

63.0 

1.306 

0.72 

1.09 

2.3 

2.5 

61.6 

2.306 

0,73 

1.02 

4.1 

1.7 

60.1 

4.11 

0.74 

0.99 

7.4 

1.2 

59.6 

7.40 

1.08 

*\ 

1.39 

* 

Adjusted  values,  see  text. 


The  median  radiant  exposure  values  for  exposure  times  .37  and  .71 
might  be  considered  low  relative  to  those  adjacent  to  them  as  shown  in 
Figure  1.  This  could  suggest  that:  (1)  for  these  two  exposure  times 
the  initial  skin  temperature  of  the  animals  was  higher  than  the  other 
animals,  thereby  requiring  less  total  energy  to  raise  their  skin  tem¬ 
perature  throng,,  the  appropriate  range  and/or  (2)  errors  in  dosimetry 
yielded  lower  measured  values  of  irradiance  than  actually  were  used. 
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A  judgment  as  to  whether  these  two  points  are  out  of  line  can  be 
made  by  considering  the  calculated  temperature  responses.  A  priori, 
it  can  be  expected  in  a  system,  such  as  the  one  considered  here,  that 
for  a  fixed  damage  response,  as  the  exposure  time  is  decreaeed,  the 
maximum  temperature  must  increase  in  order  to  compensate  for  the  short¬ 
er  time.  A  plot  of  the  maximum  temperature  values  vers  s  the  exposure 
time  should  reveal  a  smooth,  continuous  function.  The  calculated  maxi¬ 
mum  temperature  (column  3,  Table  1)  versus  the  exposure  time  is  shown 
in  the  semilogarithmic  plot  of  Figure  2.  It  can  be  seen  that  all  the 
points  cluster  closely  around  a  single  line  with  the  exception  of  the 
two  in  question.  A  calculated  18%  increase  in  irradiance  values  was 
necessary  to  give  the  maximum  temperature  values  for  these  two  points 
shown  in  Figure  1  as  open  circles.  Such  a  difference  attributed  to 
dosimetry  seems  rather  high.  The  adjusted  irradiance  values  were  used 
to  recompute  these  two  damage  integrals. 
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Fig.  2.  The  calculated  maximum  temperature  attained 
at  the  base  of  the  epidermal  level  as  a  function  of 
the  exposure  time  to  CO2  laser  radiation.  The  open 
circles  are  adjusted  values  for  the  exposure  times 
.37  and  .71  sec;  see  text  for  explanation. 


The  computed  values  of  the  damage  integrals  are  presented  graphical¬ 
ly  in  Figure  3.  All  the  values  are  grouped  closely  around  the  line  with 
the  exception  of  that  for  the  longest  exposure  time.  Calculations  for 
longer  exposure  times  indicate  that  their  damage  integral  values  are 
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even  higher.  The  calculated  damage  integral  values  for  the  exposures 
of  13  and  20  sec  are  24  and  32f  respectively.  This  indicates  that  the 
model  predicts  more  severe  damage  at  these  longer  times  than  was  actually 
seen.  If  the  dermal  capillary  bed  acts  as  a  heat  sink,  during  long  ex¬ 
posures,  it  could  keep  the  temperature  below  that  predicted  by  the  mod¬ 
el.  In  many  cases  during  long  exposure  periods,  erythema  was  evident  at 
the  site  of  exposure. 


EXPOSURE  TIME  (sac) 


i  Fig.  3.  The  calculated  damage  integral  value  for 

uniform  full  thickness  epithelial  hums  as  a  func¬ 
tion  of  the  exposure  time  to  CO2  laser  radiation. 

1  There  are  two  interesting  observations  derived  from  Table  1.  As 

the  exposure  time  is  shortened,  a  greater  percentage  of  the  total  damage 
occurs  after  the  exposure  is  terminated,  and  the  maximum  temperature  is 
reached  proportionally  later  following  the  termination  of  the  exposure. 
Both  of  these  relations  can  be  influenced  by  the  thermal  constants  of 
the  tissue  under  consideration  and  the  depth  at  which  the  temperature 
is  calculated. 

It  should  be  emphasized  that  the  numerical  values  derived  for  the 
inactivation  rate  constants  are  not  necessarily  the  true  values.  Their 
derivation  is  cased  on  the  assumption  that  the  heat  flow  equations  used 
in  calculating  the  temperature  response  in  the  tissue  accurately  describe 
the  real  situat  on.  Only  accurate  inactivation  rate  constants  would  be 
valid  for  determining  threshold  levels  for  other  situations  where  the 
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energy  absorption  pattern  is  different,  such  as  when  radiant  energy  of 
different  wavelengths  is  used  and  in  pigmented  skin. 

CONCLUSIONS 


An  analysis  of  experimental  dat3  defining  the  threshold  values  of 
irradiance  as  a  function  of  exposure  time  to  CO2  laser  radiation  to  in¬ 
duce  uniform  full  thickness  epithelial  burns  demonstrated  that  within 
the  exposure  time  range  of  0.2  to  5  sec  the  data  are  accurately  de¬ 
scribed  by  the  damage  integral  model  and  can  be  quantitatively  deter¬ 
mined  by  the  equation 


n  =  1.17  x  10 


64 


e- 49 ,11 0/ CTt 


273) 


dt 


wnere  Tt  is  the  temperature  at  the  epidermal-dermal  junction  at  time  t 
during  the  thermal  episode.  For  times  in  excess  of  5  sec  the  values  pre¬ 
dicted  by  the  model  were  higher  than  those  determined  experimentally. 

The  validity  of  the  time- temperature  relations  in  the  epidermis 
predicted  by  the  opaque  semi-infinite  thermal  model  and  the  therma1 
constants  for  porcine  skin  need  to  be  determined  and  the  model  and  con¬ 
stants  modified  if  necessary. 

Additional  experiments  to  determine  damage  thresholds  at  various 
depths  are  needed  to  extend  the  use  of  the  model  in  predicting  more  lev¬ 
els  of  damage.  Such  experiments  would  also  permit  further  evaluation  of 
some  of  the  assumptions  inherent  in  the  model. 
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